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Abstract. We used the Wide Field Planetary Camera 2 
aboard the Hubble Space Telescope to search for globular 
clusters in the inner regions of the nearby giant ellipti- 
cal galaxy NGC 5128. This galaxy is believed to be the 
product of a merger between a large elliptical galaxy and 
a small late-type spiral between 160 and 500 Myr ago. 
We identified 21 globular cluster candidates and measured 
their core radii, tidal radii, half-mass radii, ellipticities, po- 
sition angles, and V — I colors. We find evidence that the 
NGC 5128 globular cluster candidates are systematically 
more elliptical than are those of the Milky Way. Approxi- 
mately half of the candidates have {V — I)^ colors that are 
consistent with their being either old, unreddened globu- 
lar clusters, similar to those found in the Milky Way, or 
young, reddened globular clusters that may have formed 
during the recent merger event. Most of the rest have col- 
ors that are consistent with their being old globular clus- 
ters similar to those found in the Milky Way. Wc find one 
blue object with (V-I)^ < 0.26 ± 0.09. The color, red- 
dening, and integrated magnitude of this object are con- 
sistent with its being a small globular cluster with an age 
of ^ 100 Myr and a mass (based on its integrated lumi- 
nosity) of < 4000 Mq . We find no evidence for bimodality 
in the colors of the globular cluster candidates in our sam- 
ple beyond what can be explained by uncertainties in the 
differential reddening. 

Key words: galaxies: individual (NGC 5128) - galaxies: 
peculiar - galaxies: star clusters - methods: data analysis 



1. Introduction 

Globular star clusters (GCs) are among the oldest stellar 
systems in the Universe and provide a powerful tracer of 
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the evolutionary history of a galaxy. There is strong evi- 
dence that massive star clust ers ca n form during g alacti c 
mergers (e.g. Zepf & Ashman 1993 ; Schwcizcr et al. 1996 ), 
so galaxies that have recently experienced a merger event 
are ideal places to search for young GCs. Recently sev- 
eral candidates for luminous young GCs have been iden- 
tified in va rious merging galaxies, suc h as NGC 3597 
(Lutz |1991D , NGC 1275 (H oltzm an et al. |l992| ), and NGC 



7252 (Schweizer & Seitzer |l993D 



NGC 5128 {— Centaurus A, see Israel |1998| for a recent 
review) is classified as a giant SOpec galaxy. It is composed 
of a large, dominant spheroid, which itself resembles an 
EO galaxy, and a disk that contains large amounts of gas 
and dust. Soria et al. ( 1996 ) used direct observations of re- 
solved stars at the tip of the red giant branch in NGC 5128 
to determine a true distance modulus of /Xq = 27.8 ± 0.2, 
which corresponds to 3.6±0.2 Mpc, making NGC 5128 the 
nearest giant elliptical galaxy to our own. There is strong 
evidence (see the review by Ebneter & Balick 1983 ) that it 
is the product of a recent merger between a large elliptical 
galaxy and a small late-type spiral. A thick dust band is 
seen across the center of NGC 5128 and there is evidence 
for significant star formation occurring in the central re- 



gions of the galaxy. G. Harris et al. (1999) used Hubble 
Space Telescope (HST) Wide Field Planetary Camera 2 
(WFPC2) images to obtain a color-magnitude diagram 
for the outer halo of NGC 5128. They found a distance of 



3.9 ±0.3 Mpc, consistent with the Soria et al. (1996) esti- 
mate, and a population of old stars with iron abundances 

between [Fe/H] 2 and [Fe/H] ~ -hO.2. Their metal- 

licity distribution function is consistent with two bursts 
of star formation. The first having [m/H] = —0.6 and 
producing approximately one-third of the stars, and the 
second having [m/H] = and producing approximately 
two-thirds of the stars. They argue that the second burst 
of star formation occurred at least 1-2 Gyr after the first. 

The first observation of a GC in NGC 5128 was by 



Graham & Philips (1980). The galaxy i s now believed 
to have - 1700 GCs (H. Harris et al. |l984|) , wi th 87 



confirmed spectroscopically (see H. Harris et al. [1988 
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Sharpies |1988|). Recently G. Harris et al. ( |1998D used HST 
WFPC2 images to construct a color-magnitude diagram 
for C44, a GC in the halo of NGC 5128. They found that 
this GC was an old, intermediate-metallicity object simi- 
lar to the GCs in the Milky Way. G. Harris et al. ( |1992| , 
hereafter referred to as HG92) used Washington CMT1T2 
photometry to derive metallicities for 62 of confirmed 
GCs in NGC 5128 and found a mean iron abundance of 
[Fe/H] = -0.8 ± 0.2, which suggest that the NGC 5128 
GC system is ^ 3 times more metal rich than the Milky 
Way GC system. They found no evidence for any GCs 
having metallicities significantly greater than those found 
in the Milky Way GCs. Such metal-rich GCs might be ex- 
pected if some of the NGC 5128 GCs had formed recently 
in a gas-rich merger event. HG92 do, however, suggest 
that several blue GCs in NGC 5128 may be analogues of 
the intermediate-age GCs found in the Magellanic Clouds. 
On the other hand, Zepf & Ashman ( 1993 ) suggest that 
the metallicity distribution of the NGC 5128 GCs is bi- 
modal, with the high-metallicity peak at [Fe/H] = -1-0.25 
due to GCs formed in a merger. Hui et al. ( 1995| ) analyzed 
the kinematics of the NGC 5128 GC system and found 
that the metal-rich GCs are part of a dynamically sepa- 
rate system from the metal-poor GCs. Numerical simula- 
tions suggest that the merger event occurred between 160 



Table 1. Log of the observations. 



(Quillen et al. |1993| ) and 500 {D/b Mpc) Myr ago, where 
D is the distance to NGC 5128 in Mpc (Tubbs |1980D . 
This suggests that any GCs that formed in this particular 
merger should be quite young and, therefore, rather blue 
(~ 0.4 < V-I < 0.6; see Sect. |). 

Minniti et al. ( |1996| ) and Alonso & Minniti ( |1997 , 
hereafter referred to as AM97) used HST Wide- 
Field/Planetary-Camera 1 (WF/PC-1) images, taken be- 
fore the corrective optics package was installed in 1993, 
to search for GCs in the inner regions of NGC 5128. 
They identified 125 GC candidates, young associations, 
and open cluster candidates in the inner three kpc of 
NGC 5128. They also used ground-based RK photometry 
to estimate metallicities for 47 GC candidates. Schreier 



et al. ( 1996 ) found 74 compact sources along the northern 
edge of the NGC 5128 dust lane using HST WF/PC-1 im- 
ages. They estimate that most of these sources are young 
stars (spectral class A or earlier) but note that some are 
resolved and may be GCs. 

Identifying GC candidates in the inner regions of NGC 
5128 is difficult since there is nonuniform extinction, con- 
tamination from foreground stars and background galax- 
ies, and confusion with open clusters and blue, star- 
forming knots in NGC 5128. GC candidates can not 
be identified based solely on their colors since the large 
amount of uneven reddening makes it very difficult to 
determine the dereddened color of an object. A better 
approach is a scheme to identify GC candidates based 
solely on their structural parameters. All known GCs in 
the Local Group can be r eason ably w ell fit by Michie- 
King models (Michie |l963| ; King |l966|) , ahhough - 20% 
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show evidence of having undergone core collapse. The vast 
majority of the Milky Way's GCs are uniformly old ob - 
jects with ages of 11.5 ± 1.3 Gyr (Chaboyer et al. 1998 ), 
mean King core radii of re — 2.3 ± 0.4 pc, mean King 
tidal radii of rt — 45.2 ± 3.5 pc, mean concentrations of 



c = logiQ{rt/rc) = 1.40 ± 0.04, and mean elhpticities of 
e = 0.08 ± 0.01. If the GCs in NGC 5128 are structurally 
similar to those in the Local Group spiral and dwarf galax- 
ies, then high resolution imaging can be used to identify 
GC candidates in the inner regions of NGC 5128 with- 
out resorting to an identification scheme based upon the 
integrated colors of the objects. 



2. The Data 

2.1. Observations 

We used the WFPC2 aboard the HST to obtain F555W 
(WFPC2 broadband V) and F814W (WFPC2 broadband 
/) images of a region near the nucleus of NGC 5128. The 
data were obtained on July 27, 1997 for the cycle 6 pro- 
gram GO-6789. The WFPC2 had an operating tempera- 
ture of — 88°C and a nominal gain setting of 7 e~/ADU. 
The observations are listed in Table m 

Exposures were taken in each field with each of the 
F555W and the F814W filters. Cosmic rays impact ~ 20 
pixels per second on each WFPC2 CCD, but by combining 
the multiple exposures per filter for each field, the number 
of pixels lost to cosmic ray events is negligible. Therefore, 
we did not apply any processing explicitly to remove cos- 
mic rays from the images. The data were preprocessed 
through the standard STScI pipeline for WFPC2 data. 
Known bad pixels were flagged and not used in the data 
analysis. No corrections were made for geometric distor- 
tions in the area of the WFPC2 pixels. 

The survey consists of six fields that cover a to- 
tal area of approximately 5.'3 x 8f0 centered on a = 
13'^25'"27!3, ,5 = -43°01'09" (J2000 coordinates), the nu- 
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cleus of NGC 5128. Adjacent fields overlap by ~ 0.'5 giving 
a total effective area of ~ 25 D' for the survey. 

2.2. Data Reductions 

We combined the exposures for each field by taking the 
average of the three images in each filter (four images for 
the F555W exposures of Field 2). No re-registration of the 
images was performed since the shifts between the images 
were typically less than 0.1 pixel (O'.'Ol on the WFC and 
0'.'005 on the PC). We estimate that combining the images 
in this way may result in the sizes of the GC candidates be- 
ing systematically overestimated by no more than ~ 0'.'02. 
We prefer to introduce this simple systematic offset than 
deal with the poorly-understood systematic uncertainties 
that arise from interpolating flux across fractional-pixel 
shifts. 

2.2.1. Identifying Globular Cluster Candidates 

At the distance of NGC 5128 (d = 3.6±0.2 Mpc), the mean 
King core- and tidal-radii of the Milky Way GCs would 
appear to be f^ = 0'.'13±0'.'02 and n = 2'.'59±0'.'20, respec- 
tively. Therefore, any GCs in NGC 5128 will appear to be 
semi-stellar and be strongly affected by the point spread 
function (PSF) of the WFPC2. After some experimenta- 
tion, we adopted the following procedure for identifying 
GC candidates. We wish to stress that this procedure is 
quite strict and will probably result in the rejection of 
some legitimate GC candidates. However, we prefer to re- 
ject real GCs rather than have our sample contaminated 
with stars or background galaxies. 

In order to increase the signal-to-noise ratio (S/N) of 
the GC candidates - a particularly important point for the 
faint {V < 20) GC candidates - we combined the F555W 
and F814W images for each field to get finding images. 
The dust lane introduces variations in the background on 
spatial scales of ~ 1", comparable to the expected sizes 
of the GC candidates in NGC 5128. To reduce the effects 
of the uneven background light, large-scale spatial varia- 
tions in the background were removed by running a ring 
median filter (Seeker 1995) over the finding image, sub- 
tracting the resulting smoothed background, and adding 
back the mean background value. The median filter radius 
was set to 1", which is ^ 3.5 times the expected full- width 
at half maximum (FWHM) of a typical GC candidate. 
This choice of filter radius ensures that the cores of the 
GC candidates will not be altered by the median filter and 
that any background structure larger than a typical GC 
candidate will be removed. Since the most extended Milky 
Way GCs have tidal radii that are significantly greater 
than 2.5 times their FWHM, and extended halos have 
been detected around several Galactic and extra-Galactic 



GCs (GriUmair et al. 1995; 1996; Holland et al. 1997) 



this approach will alter the distribution of light in the 
outer regions of most of the GC candidates. However, this 



is not important since the finding images are used only 
to construct a preliminary list of GC candidates. A more 
rigorous set of criteria, based on the structures of the GC 
candidates as determined from the original images, will 
be applied to the preliminary list to obtain a final list of 
GC candidates in the central regions of NGC 5128. 

The first step in our identification procedure was to 
run the DAOPHOT ii (Stetson 1987 ; 1994 ) find routine 
on the background-subtracted images to identify GC can- 
didates. The finding thresholds were set to Qasky for the 
PC images and lOtTsky for the WFC images. Tests with 
artificial GCs suggested that any detections below these 
thresholds would be rejected at some point in our identifi- 
cation process. Daophot ii find has an algorithm for re- 
jecting non-stellar objects based on two parameters called 
"sharpness" and "round" . This algorithm was turned off 
since images of GCs can have different shapes and concen- 
trations from images of stars. 

Next, the daophot ii photometry routine was used 
to obtain aperture photometry for each of these detec- 
tions. The photometry was performed separately on each 
of the combined F555W and F814W frames, not on the 
combined finding frame. An aperture radius of 0'.'2 was 
used since most Galactic GCs, if moved to the distance 
of NGC 5128, would appear to have core radii smaller 
than this. Therefore, the signal within the aperture will be 
dominated by the light from the object and not from the 
background. Candidate objects with S/N < 5 within the 
photometry aperture were discarded since the signal was 
not strong enough to determine reliable shape parameters. 
The sky brightness was determined in an annulus with an 
inner radius of ~ 0'.'9 and an outer radius of ~ I'.'l. This 
annulus was chosen to be far enough from the center of 
the GC candidate that the light in the annulus will be 
dominated by the background, yet near enough to the GC 
candidate that the light in the annulus will be a reasonable 
approximation of the mean background at the location of 
the object. For large GC candidates this annulus will be 
inside the tidal radius of the object so our estimate of 
the background will be contaminated. However, the values 
determined at this stage are only preliminary estimates, 
which will be improved upon later in the identification 
process when Michie-King models are fit to the GC can- 
didates. The lists of GC candidates in each of the F555W 
and F814W images were matched using the daomatch 
and DAOMASTER software. Only objects that appeared in 
both the F555W and F814W images, and whose centers 
matched to within 0'.'05 {^ 1.1 pixel on the PC images 
and '~ 0.5 pix on the WFC images), were considered to be 
real GC candidates. 

Distinguishing bona fide GCs from stars and back- 
ground galaxies is challenging. The colors of the objects 
can not be used since we are interested in studying the 
color distribution of GCs in NGC 5128 and do not wish to 
bias our sample. To make matters worse, the presence of 
dust in NGC 5128 will add a significant amount of scatter 
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to the intrinsic color distribution, and may cause legit- 
imate GCs to be rejected if a color-based identification 
scheme is used. The solution is to identify GC candidates 
by their structural parameters, although the best choice 
of structural parameters is not obvious. At the distance 
of NGC 5128 a typical Galactic GC would appear to have 
an intrinsic FWHM of ~ 0'.'25, or approximately twice 
the FWHM of the WFPC2 PSFs. Therefore, the observed 
FWHM, concentration, and ellipticity of a GC candidate 
can be heavily influenced by the PSF. Since the PSF varies 
strongly with position on the WFPC2 CCDs, the poten- 
tial for confusion between stellar images and concentrated 
GC candidates is great if the PSF is not removed, in some 
way, from the data. Therefore, the observed shape of an 
object can not be directly used to classify it as a star, GC 
candidate, or galaxy. 

After some experimentation with adding and recover- 
ing artificial GCs and artificial stars, we found that the 
following procedure was reasonably reliable for identify- 
ing GC candidates. For each GC candidate we took all 
the pixel values within 1" of the center of the object and 
subtracted an estimate of the local background (the center 
and background were determined by the daophot ii pho- 
tometry algorithm) . A one-dimensional Moffatian (Mof- 
fat |l969l). 



Mir,s) = M(0) 



' + (v) 



(1) 



was fit to each candidate using the effective radius, roff , in- 
stead of the true distance from the center of the candidate 
in order to compensate for any ellipticity that might be 
introduced by the PSF. The effective radius of an ellipse 
is defined by abir = r^fjTr, where a and b are the lengths of 
the semi-major and semi-minor axes of the ellipse, respec- 
tively. For each pixel the effective radius from the center 
of the GC candidate was computed using: 



'eff 



(1 - e) cos^ 6*0 + sin^ 6lo /(I - e) 



cos^ 6lo -f- (1 - e) sin^ 6*0 /(I - e) 



-|- 2xye{e — 2) cos6'osin0o/(l ~ e), 



(2) 



where x and y are the coordinates of the pixel on the CCD, 
and e and ^o a-re estimates of the ellipticity and position 
angle of the GC candidate. The latter two quantities were 
estimated from the moments of the light from each object. 
In order to determine which combinations of a and 
P corresponded to stars and which corresponded to GC 
candidates, a series of artificial stars and artificial GCs 
were added to the images. The artificial stars were added 
using the DAOPHOT n ADDSTAR routine and the appropri- 
ate PSFs scaled to magnitudes of 16 < F555W < 22. The 
artificial GCs, also with integrated magnitudes of 16 < 
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Fig. 1. This figure shows the best-fitting Moffatian a and 
P parameters for the artificial stars (circles) and artificial 
GCs (crosses). Based on the distribution of objects in this 
diagram we assumed that any objects that lie inside the 
wedge formed by the solid lines were GC candidates. No 
upper limit was placed on the value of /?. 



F555W < 22, were added using the IRAF^ v2.10.4 task 
NOAO.ARTDATA.MKOBJECT. The artificial GCs all had el- 
lipticities of e = 0, concentrations of 0.67 < c < 2.12, and 
core radii of 0'.'067 < Vc < 0'.'4 (corresponding to physi- 
cal core radii of between ~ 1 and ~ 7 pc at the distance 
of NGC 5128). Therefore, the artificial GCs had a range 
of structures similar to those of the Milky Way's GCs. 
The procedure described above was used to determine the 
Moffat a and (3 for each artificial object. The results are 
presented in Fig. m and were used to determine which com- 
binations of a and /3 represent stars and which represent 
GC candidates. These limits on a and f3 were then ap- 
plied to the GC candidates found on the WFPC2 images 
of NGC 5128. 

Figs. and show the Moffatian a and P parame- 
ters for the ~ 3800 objects that were successfully fit by 
Moffatians. There are 403 objects with Moffat parameters 
that lie inside the wedge (see Fig. p. Our simulated data 
suggest that these are extended objects such as GCs, back- 
ground galaxies, dust features, open clusters, star forming 
regions, or blended stars. 



^ Image Reduction and Analysis Facility (IRAF) , a software 
system distributed by the National Optical Astronomy Obser- 
vatories (NOAO). 
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Fig. 2. This figure shows the best- fitting Moffatian a and 
P parameters for each object (small circles) on the F555W 
images. Simulated data (see Fig. |^) suggest that objects 
that lie inside the wedge formed by the solid lines are ex- 
tended objects. Therefore we consider any objects that lie 
inside the wedge to to be GC candidates. The solid squares 
show the locations of the GC candidates from Table 0. 
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Fig. 3. This figure shows the best-fitting Moffatian a and 
P parameters for each object (small circles) on the F814W 
images. Simulated data (see Fig. |l|) suggest that objects 
that lie inside the wedge formed by the solid lines are ex- 
tended objects. Therefore we consider any objects that lie 
inside the wedge to to be GC candidates. The solid squares 
show the locations of the GC candidates from Tabic |[ 



2.2.2. Fitting Michie-King Models 

We fit a two-dimensional, PSF-convolved, single-mass 
Michie-King model to each of the 4 03 G C candidate us- 
ing software developed by Holland ( 1997 ). This software 
assumes that the surface brightness profile along the ef- 
fective radius axis of a GC candidate with an ellipticity 
of e and a position angle of 9o has a King profile with a 
concentration of c and a core radius of Tc- It then builds 
a two-dimensional model based on this surface brightness 
profile, e, and ^o- The two-dimensional model is convolved 
with the appropriate PSF for the location on the CCD and 
a chi-square minimization is performed between the PSF- 
convolved model and the original data image. The soft- 
ware uses CERN's minuit function minimization package 
to fit simultaneously the concentration, core radius, to- 
tal flux in the object, ellipticity, position angle, and mean 
background. Objects located within 32 pixels of the edge 
of a CCD (= 3'.'2 for the WFC and 1'.'6 for the PC) were 
not fit to avoid the edges of the CCD biasing the fits. Once 
a best fit had been determined, the King tidal radius, rt, 
and the half-mass radius, rh, of the model were computed. 
Separate fits were made to the F555W images and the 
F814W images and an object was considered to be GC 



candidate only if a Michie-King model could be fit in both 
colors. We were able to fit Michie-King models to 98 of the 
403 potential GC candidates. Mean structural parameters 
were calculated for these object by taking the mean of 
the values found in each filter. Four objects (#8, #113, 
#128, and #129) (see Tables |, |, and |) were identified 
on multiple fields. In these cases we computed the mean 
of the structural parameters measured in each field. 

We elected to separate GC candidates from back- 
ground galaxies based on their fitted ellipticities and half- 
mass radii (see Fig. ^). Half- mass radii are preferred to 
tidal radii or core radii because Fokker-Planck models 
of spherical stellar systems show that half-mass radii re- 
main reasonably constant over periods of several Gyr 
(e.g. Cohn 1979; Takahashi 1997), making it a unique 
length scale for GCs. The mass interior to the half- 
mass radius tends to undergo a gravo-thermal collapse 
and become concentrated at the center of the GC over 
time (i.e. core-collapse), which results in the core radius 
shrinking. Meanwhile, the mass exterior to the half-mass 
radius tends to expand outwards, causing the tidal ra- 
dius to grow. Since we are interested in finding young, 
intermediate-age, and old GCs in NGC 5128, it is useful to 
have a selection criterion that does not depend on the age 
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Fig. 4. The ellipticity vs. half-mass radius of the best- 
fitting single-mass Michie-King model for each object 
where a Michie-King model was successfully fit. Objects 
with rh > 10" {^ 175 pc) have not been plotted. The solid 
box in the lower left of the plot shows the region occupied 
by Galactic GCs. Based oir this plot we have assumed that 
any object with rh < 2" {^ 35 pc) and e < 0.4 (the dashed 
box) is a GC candidate in the NGC 5128 system. 



of the GC candidate. Galactic GCs have half-mass radii of 



approximately 1.3 < r^ < 31.9 pc (W. Harris 1996 ), which 
corresponds to 0'.'07 < Vh < 1'.'83 at the distance of NGC 
5128. There is no evidence that the radius of a Galactic GC 
depends on its mass (van den Bergh et al. 1991). There- 



fore, we have assumed that only objects with rh < 2" 
(~ 35 pc at the distance of NGC 5128) were GC candi- 
dates. It is possible that some of the objects in Fig. |^ that 
have high ellipticities and low half-mass radii are double 
clusters. However, Innanen et al. 1983 have shown that 
a binary GC could not survive a single Galactic orbit in 
the Milky Way so it is unlikely that there are any old, 
or intermediate-age double GCs in NGC 5128. It is possi- 
ble that very young multiple GCs that formed within the 
last ^ 100 — 200 Gyr could have survived to the present 
day, but we are unable to differentiate between them and 
background galaxies. 

The most elliptical Galactic GC is M19 with e = 
0.27 (White & Shawl p^) and the most elliptical GC 



known is NGC 2193 in the Large Magell anic C loud (LMC) 
which has e = 0.33 (Geisler & Hodge |1980D. Geisler & 



Hodge (1980) modelled the distribution of observed el- 



lipticities for 25 GCs in the LMC and found that it was 
unlikely that the largest true ellipticity exceeded e = 0.4. 
The LMC contains both dynamically young and dynami- 



Fig. 5. A finding chart for Field 1. The GC candidates 
are circled with their identification numbers (see Table ^ 
printed near each object. 

Fig. 6. A finding chart for Field 2. 

Fig. 7. A finding chart for Field 3. 

Fig. 8. A finding chart for Field 4. 



cally old GCs, so the largest ellipticity seen in the LMC is 
a reasonable estimate of the largest ellipticity that we can 
expect to see in NGC 5128. Therefore, only objects with 
e < 0.4 were considered to be GC candidates. 

The final step was to examine visually the WFPC2 
images of each GC candidate to ensure that the Michie- 
King model fits looked realistic. We found that ~ 20% 
of the objects were either located on diffraction spikes 
from saturated stars, or exhibited unusually large resid- 
uals when the best-fitting Michie-King models were sub- 
tracted. These spurious identifications were discarded. 

Fig. ^ shows the measured half-mass radii and ellip- 
ticities for the surviving objects in NGC 5128 and Ta- 
ble H shows the final list of GC candidates that we find in 
the central regions of NGC 5128. The second and third 
columns show the J2000 coordinates of the objects as 
determined using the IRAF/STSDAS (v2.0.1) task STS- 
DAS. TOOLBOX. imgtools.xy2rd. Column 4 is the ob- 
served (projected) distance of the GC candidate from the 
center of NGC 5128 in arcminutes. The center of NGC 
5128 was taken to be aj2ooo = 13h25'"27;3, (5j2ooo = 



-43°01'09" (Johnston et al. |1995| ). Columns 5 and 6 give 
the field (from Table m and CCD that the object was 
found on. Columns 7 and 8 give the X and Y coordinates 
(in pixels) on the CCD. Column 9 lists the identification 
number of the object in Table 1 of Minniti et al. (1996). 

Tables y and U lists the coordinates for the 61 extended 
objects with rh > 2' and e > 0.4. Some of these objects 
may be GCs in NGC 5128 while others may be background 
galaxies with structures similar to those of Michie-King 
models. Six of these objects have been previously identi- 
fied as GCs by Minniti et al. (|l996| ) and Sharpies ( |1988[ ). 

Figs. pI through pi]| show the locations of the 21 GC can- 
didates on the F814W-band WFPC2 images. Only objects 
that pass all of the criteria described above are marked on 
these figures. Objects (such as #15) were only marked on 
the fields that they were identified as GC candidates in. 
In most of the cases where a GC candidate is present in 
multiple fields, but only identified in one field, the GC 
candidate was located very near the edge of one of the 
CCDs. Spatial variations in the PSF are largest near the 
edges of the CCDs so the Michie-King model fits are less 
reliable near the edges of the CCDs. 
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Table 2. The GC candidates in the central regions of NGC 5128. 



ID 



QJ2000 



D Field CCD 



X 



Y 



Other 



1 


13'^25'"16'.'1 


-42°59'45" 


2.'50 


3 


3 


703.990 


570.660 




2 


13''25"16'.'4 


-43°02'10" 


2. '26 


6 


3 


726.813 


361.345 




3 


13''25'"18'.'1 


-43°01'59" 


1.'90 


6 


3 


516.154 


406.027 




4 


13''25'"19'.'4 


-43°01'13" 


1.'47 


3 


2 


482.000 


599.000 




5 


13''25'"20'.'5 


-43°00'52" 


1.'30 


3 


2 


310.789 


429.262 




6 


13''25'"25'.'3 


-43°02'01" 


0.'95 


5 


3 


341.557 


317.009 




7 


13''25'"26'.'6 


-42°59'27" 


1.'71 


1 


3 


686.555 


752.859 


M27 


8 


13''25'"27'.'0 


-43°00'02" 


1.'12 


1 


3 


746.948 


397.874 


M26 


8 


13''25"'27'.'0 


-43°00'02" 


1.'12 


2 


3 


171.527 


228.000 


M26 


9 


13''25'"27'.'1 


-42°59'40" 


1.'49 


2 


3 


101.013 


438.000 


M25 


10 


13''25™27'.'2 


-42°59'28" 


1.'69 


1 


3 


621.989 


725.544 


M24 


11 


13''25'"27'.'2 


-43°01'54" 


0.'75 


4 


3 


672.131 


515.204 




12 


13''25'"27'.'8 


-42°59'23" 


1.'77 


1 


3 


551.791 


755.891 


M23 


13 


13''25""28'.'5 


-43°02'57" 


1.'81 


5 


2 


432.822 


211.852 




14 


13''25'"29'.'2 


-43°00'16" 


0.'94 


1 


3 


553.540 


198.363 


M20 


15 


13''25'"29'.'8 


-43°00'07" 


1.'12 


1 


3 


463.865 


268.523 


M18 


16 


13''25'"30'.'2 


-42°59'36" 


1.'63 


1 


3 


337.483 


548.311 


M17 


17 


13''25'"31'.'1 


-42°59'39" 


1.'65 


1 


3 


243.088 


494.704 




18 


13''25'"31'.'5 


-43°00'04" 


1.'32 


1 


3 


280.111 


236.550 


M15 


19 


13''25'"31'.'7 


-43°00'35" 


0.'97 


1 


2 


160.847 


343.070 




20 


13''25"33'.'7 


-43°01'20" 


1.'16 


1 


2 


655.839 


265.242 




21 


13''25'"34'.'4 


-43°03'30" 


2.'67 


4 


2 


763.000 


276.549 





Fig. 9. A finding chart for Field 5. 
Fig. 10. A finding chart for Field 6. 

The structural parameters of the best-fitting Michie- 
King models, as well as the fitted ellipticities and position 
angles, for each GC candidate are listed in Table 0. The 
first column contains the object ID (from Table ^. The 
various radii are given in units of seconds of arc, and the 
position angles are measured in degrees with 9o = 0° being 
north and 6o increasing to the east. The xl values are the 
reduced goodness of fit values returned by the fitting soft- 
ware. The uncertainties (a) are the standard deviations in 
the values for the parameters that were measured in each 
filter. The position angles for GC candidates with small 
ellipticities (< 0.05) are not reliable. All of the xt are sig- 
nificantly less than one, which suggests that the formal 
uncertainties in the model's parameters are not reliable. 
Therefore, we have elected to estimate the uncertainties 
in the fits through monte-carlo simulations as described 



in Sect. 2.4 



2.3. Contamination 

Our data will contain images of Galactic foreground stars, 
and supe rgian ts in NGG 5128. From the work of Bahcall & 
Soneira ( 1981 ) we expect to find ~ 700 Galactic stars with 
/ < 20 in our fields. The brightest stars in the halo of NGC 



5128 have V - 24.5 (Soria et al. |1996D while the bright- 
est young stars in the central regions of NGC 5128 are 
expected to have approximately 18 < V^ < 21. The mor- 
phological selection criteria that we applied to our data 
are very effective at rejecting stars (see Figs. |], ||, and ^,), 
and the HST images show all of the 21 GC candidates to 
be extended objects, so we believe that stellar contamina- 
tion is not a problem in our data. 



The galaxy counts of Tyson (1988) suggest that there 
will be '^ 10 background galaxies in our images down to 
/ ^ 20. However, many of these galaxies will be obscured 
by the dust lane, so we will detect significantly fewer than 
this. The morphological criteria that we applied to obtain 
our list of GC candidates will reject any galaxies that are 
not morphologically similar to the GCs found in the Milky 
Way or LMC. AM97 used a comparison field located 30' 
northeast of the nucleus to estimate that ^ 20% of the ob- 
jects that they detect in their search for GCs in the inner 
regions of NGC 5128 are foreground stars or background 
galaxies. Since the morphological criteria that we applied 
are stricter than those of AM97, we believe that 20% is a 
reasonable upper limit on the amount of contamination in 
our list of GC candidates. 



2.4- Uncertainties in the Structural Parameters 

The MINUIT package provides an estimate of the formal 
uncertainty in each parameter based on the covariance 
matrix of the fit. In general the formal uncertainties were 
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Table 3. Extended objects with rh > 2" and e > 0.4 in the central regions of NGC 5128 



ID 


Q.J2000 


Sj2000 


D 


Field 


CCD 


X 


Y Othei 




101 


13*^25'" 19'.' 3 


-43°02'30' 


2:00 


6 


3 


508.484 


64.127 G331 


102 


13''25'"20'.'0 


-43°01'56' 


1:57 


6 


3 


307.668 


352.171 


103 


13''25'"20'.'6 


-43°01'11' 


1:25 


3 


2 


500.453 


470.957 




104 


13''25'"21'.'0 


-42°59'25' 


2:10 


3 


3 


128.163 


607.910 




105 


13''25'"21'.'8 


-43°01'26' 


i;07 


6 


4 


561.226 


61.192 




106 


13''25'"22'.'1 


-43° 00' 39' 


1:10 


3 


2 


234.000 


214.221 




107 


13''25'"22'.'7 


-43° 00' 39' 


l!00 


3 


2 


256.978 


156.827 




108 


13''25'"22'.'8 


-43°00'40' 


0:98 


2 


2 


99.928 


720.877 




109 


13''25'"23'.'0 


-43°01'18' 


0!82 


6 


4 


593.184 


215.648 




110 


13''25'"23'.'1 


-43°00'50' 


o:85 


3 


2 


377.582 


144.000 




111 


13''25'"23'.'3 


-43°01'48' 


0:99 


5 


3 


497.984 


517.934 




112 


13''25'"23'.'6 


-43°01'34' 


o:8i 


5 


3 


421.420 


632.818 




113 


13''25'"23'.'7 


-43°00'47' 


0^78 


2 


2 


195.364 


642.153 




113 


13''25'"23'.'7 


-43°00'47' 


0:78 


3 


2 


369.000 


69.211 




114 


13'^25'"24'.'5 


-43°01'05' 


0:54 


6 


4 


659.413 


414.000 




115 


13''25'"24'.'9 


-43°00'20' 


0:94 


2 


3 


444.626 


122.783 




116 


13''25'"25'.'0 


-42°59'40' 


1:55 


2 


3 


326.147 


506.868 




117 


13''25'"25'.'3 


-43°02'01' 


0:95 


6 


4 


101.986 


296.981 




118 


13''25'"25'.'6 


-43°00'58' 


0:38 


2 


2 


354.000 


478.498 




119 


13''25'"25'.'9 


-43°00'57' 


0f35 


2 


2 


358.270 


438.501 




120 


13''25'"26'.'0 


-43°02'16' 


1:14 


6 


1 


629.884 


296.782 




121 


13''25"'26'.'5 


-42°59'55' 


1:25 


2 


3 


207.559 


311.875 




122 


13''25'"26'.'6 


-42°59'27' 


1:71 


3 


4 


412.549 


539.948 M27 




123 


13''25'"26'.'6 


-43°02'11' 


i:o4 


5 


3 


241.413 


170.644 




124 


13''25'"26'.'8 


-43°00'12' 


0:97 


2 


3 


224.417 


144.083 




125 


13''25'"26'.'9 


-43°00'14' 


0:93 


2 


3 


214.154 


115.666 




126 


13''25'"27'.'1 


-43°00'06' 


i:o6 


1 


3 


744.088 


355.292 




127 


13''25'"27'.'1 


-43°01'54' 


0f75 


5 


3 


125.479 


305.598 




128 


13''25'"27'.'3 


-42°59'28' 


1:59 


2 


3 


45.380 


554.534 M24 


128 


13''25'"27'.'3 


-42°59'28' 


1:59 


3 


4 


386.949 


604.233 M24 


129 


13''25'"27'.'3 


-43°01'54' 


0:75 


5 


3 


113.000 


300.561 


129 


13''25'"27'.'3 


-43°01'54' 


0:75 


6 


4 


87.553 


525.160 




130 


13'^25'"27'.'4 


-43°01'28' 


0:31 


2 


2 


711.369 


371.000 




131 


13''25'"27'.'7 


-43°00'13' 


0:93 


1 


3 


698.887 


273.893 




132 


13''25'"27'.'8 


-42°59'23' 


lf77 


2 


4 


582.795 


98.323 




133 


13''25'"28'.'3 


-43°00'19' 


o:85 


1 


3 


658.173 


199.939 




134 


13''25'"28'.'3 


-43°01'38' 


0:51 


5 


4 


408.000 


133.199 




135 


13''25'"28'.'5 


-43°01'23' 


0:30 


5 


4 


536.658 


206.354 




136 


13''25'"29'.'0 


-43°00'43' 


0:52 


2 


2 


323.925 


78.964 




137 


13''25'"29'.'2 


-43°00'16' 


0:94 


2 


1 


162.995 


99.487 




138 


13''25'"29'.'2 


-43°00'22' 


o:85 


1 


3 


570.063 


142.000 




139 


13''25'"29'.'2 


-43°01'28' 


0^45 


4 


3 


374.195 


675.000 




140 


13''25'"29'.'2 


-43°02'29' 


1:37 


4 


3 


599.716 


108.340 





~ 10% of the best-fit value of each parameter. This is 
consistent with the uncertainties quoted in Table |5| that 
were determined from the differences between the best 
fitting structural parameters determined from the F555W 
images and the F814W images. 

In order to test the formal uncertainty estimates, and 
to look for systematic differences between the recovered 
structural parameters and the true structural parameters 
of the NGC 5128 GC candidates, we constructed a series 



of artificial GCs and added them to the NGC 5128 images. 
The total of 810 artificial GCs were added to the WFPC2 
images with randomly assigned concentrations between 
0.67 < c < 2.74, core radii between 0'.'07 < r^ < 0'.'4, ellip- 
ticities of e = 0, and magnitudes of 16 < F555W < 22. We 
found that the recovered concentrations for the brightest 
artificial GCs (16 < F555W < 18) were within ~ 5% of 
their true values 95% of the time. For the faintest artifi- 
cial GCs (20 < F555W < 22) the recovered values were 
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Table 4. Extended objects with r^ > 2" and e > 0.4 in the central regions of NGC 5128, continued. 



ID 


Q.J2000 


5.72000 


D 


Field 


CCD 


X 


Y 


Other 


141 


13''25'"29'.'2 


-43° 00' 38" 


o:6i 


1 


2 


109.346 


606.980 




142 


13''25'"29'.'7 


-43°01'14" 


0!42 


1 


2 


470.698 


664.338 




143 


13''25'"29'.'8 


-43°00'07" 


i:i2 


2 


4 


95.982 


190.565 


M18 


144 


13''25'"30'.'0 


-43°02'12" 


i:i5 


4 


3 


460.000 


234.000 




145 


13''25'"30'.'1 


-43°02'09" 


Ifll 


4 


3 


430.157 


264.000 




146 


13''25'"30'.'3 


-43°01'12" 


0:53 


1 


2 


466.932 


598.975 




147 


13''25'"30'.'6 


-43° 00' 24" 


0!95 


1 


3 


432.260 


75.170 




148 


13'^25'"31'.'1 


-43°01'19" 


0:59 


5 


4 


476.273 


487.763 




149 


13''25'"31'.'2 


-43°01'19" 


0!71 


1 


2 


567.825 


517.798 




150 


13''25'"31'.'3 


-43°00'37" 


o:89 


1 


2 


169.525 


384.000 




151 


13''25'"31'.'5 


-43°00'05" 


1:30 


2 


4 


66.422 


375.512 


M15 


152 


13''25'"31'.'8 


-43°01'56" 


1:12 


5 


4 


100.413 


429.480 




153 


13''25'"32'.'6 


-43°01'27" 


0f99 


5 


4 


345.530 


611.654 




154 


13''25'"33'.'0 


-42°59'05" 


2:31 


1 


4 


761.480 


131.378 


M12 


155 


13''25'"33'.'1 


-42°59'05" 


2:32 


2 


4 


595.554 


709.645 




156 


13''25'"33'.'2 


-43°01'37" 


1:15 


5 


4 


224.435 


632.127 




157 


13''25'"33'.'6 


-43°01'20" 


1:14 


5 


4 


373.122 


735.354 




158 


13''25'"34'.'6 


-43°02'37" 


1:95 


4 


2 


280.000 


74.000 




159 


13''25'"35'.'9 


-42°59'50" 


2:04 


1 


4 


233.755 


300.782 




150 


13''25'"37'.'0 


-43°01'29" 


lf78 


4 


4 


360.651 


495.008 




161 


13''25'"37'.'8 


-43°01'18" 


1:90 


4 


4 


433.330 


620.504 





within 22% of the true values 95% of the time. Systematic 
shifts were negligible for artificial GCs with concentrations 
greater than c ~ 0.9, but grew rapidly for artificial GCs 
with concentrations less than this. The uncertainties in 
the structural parameters increased as the concentration 
decreased. Other structural parameters behaved in very 
similar manners. 

A second source of uncertainty in the fitted structural 
parameters is the uncertainty in the PSF. Each GC can- 
didate was fitted with a Michie-King model that was con- 
volved by an estimate of the PSF at the location of the ob- 
ject on the CCD. We used PSFs that were created by Pe- 
ter Stetson (1996, private communication) from WFPC2 
observations of stars in the Galactic GC uj Centauri. How- 
ever, the NGC 5128 images were taken approximately 
three years after the w Gen ones, so long-term variations in 
the focus of the WFPC2 (e.g. Suchkov & Casertano |1997D 
raise the possibility of a mismatch between the actual and 
adopted PSFs. If this is the case then the derived struc- 
tural parameters would be in error, as the fitted Michie- 
King models were convolved with a PSF constructed us- 
ing a slightly different focus from that of our images. In 
order to estimate the possible errors introduced by uncer- 
tainties in the PSF, we repeated the Michie-King model 
fitting procedure with the wrong PSFs. In other words, 
we fit Michie-King models to the GC candidates on the 
WF3 CCD using both the F555W and the F814W PSFs 
from the WF2, WF3, and WF4 CCDs. This gave us six 
estimates of the structural parameters of each object ob- 
tained using six variations of the WF PSF. For each of 



these estimates we computed the difference between Wo 
as determined using the correct PSF, and the five WqS de- 
termined using the incorrect PSFs (AWq). We then com- 
puted the mean (AWq) and standard error in the mean of 
these five values for each GC candidate. This gave us an 
estimate of the systematic uncertainty in the value of Wo 
that we derived for each GC candidate. Finally, we com- 
puted the mean, standard error in the mean, and median 
of the individual AWq values for all the GC candidate. 
This gave an estimate of the systematic uncertainty for 
a typical GC candidate in our data. These values, along 
with analogous estimates of the systematic errors in the 
other structural parameters, are presented in Table g. This 
technique is not mathematically rigorous since the varia- 
tions in the PSF from one CCD to another, and from one 
filter to another, are not the same as the variations due 
to changes in the focus of the HST over a per iod o f two 
or three years. However, Suchkov & Casertano (1997) find 
that long-term changes in focus introduce changes of a 
few percent in the photometric calibrations, which corre- 
sponds to changes of a few percent in the shape of the 
PSF. The variations in the PSF from one CCD to the 
next can be much larger than this so we believe that the 
systematic errors derived here represent an over-estimate 
of the true systematic errors introduced by possible long- 
term variations in the PSF. 



3. Structural Parameters 
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Table 5. The best-fitting structural parameters for tlie NGC 5128 GC candidates. Tfie values are the means of the 
structural parameters derived from the F555W and F814W images. 



ID Wo±a 



Tc ± (7 



rh±o- 



rt±(T 



c± a 



€±(J 



6*0 ±cr 



"^ 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 



6.2±--- 
6.8 ±0.3 

4.1 ±1.1 
8.6 ±0.2 

6.2 ±4.1 

6.2 ±0.2 
3.5 ±1.2 
9.0 ±0.2 

7.3 ±0.1 

5.5 ±0.6 

6.6 ±1.0 

4.2 ±1.0 

7.3 ±0.1 
8.0 ±0.8 

4.8 ±0.1 
7.2 ±0.4 

4.7 ±1.6 

6.9 ±0.0 
5.0±--- 
7.2 ±1.0 
3.9 ±2.0 



07067 ± 
0'.'053 ± 
O'.'113±0 
0'.'024 ± 
0'.'070 ± 0' 
O'.'140±0 
O'.'149±0 
0'.'014 ± 
0'.'033 ± 

o'.'ogi ± 

O'.'102±0 
0'.'092 ± 
0'.'039 ± 
0'.'015 ± 
0'.'084 ± 0' 
0'.'044 ± 
0'.'048 ± 
0'.'062 ± 
O'.'123±0 
0'.'050 ± 
O'.'138±0 



001 
.'006 
.'035 
.'017 
.'076 
.'001 
.'101 
.'006 
.'001 
,'015 
,'001 
,'023 
,'001 
,'009 
,'006 
,'009 
,'018 
,'000 
,'010 
,'025 
,'058 



0'.'816 
^.'423 
1'.'177 
0'.'741 
1'.'428 
0'.'413 
0'.'452 
0'.'708 
0'.'667 
1'.'551 
0'.'371 
0'.'822 
0'.'460 
0'.'419 
0'.'867 
0'.'245 
0'.'982 
0'.'71 
0'.'958 
0'.'511 



± 0'.'023 
± O'.'OSO 
± 0'.'057 
± 0'.'726 
± 0'.'564 
±0'.'167 
±0'.'122 
± 0'.'078 
± 0'.'037 
±0'.'100 
± 0'.'920 
± 0'.'062 
± 0'.'029 
± 0'.'045 
± O'.'Oll 
± 0'.'041 
± 0'.'076 
± 0'.'018 
1±... 

±0'.'132 
±0'.'191 



f.'317 ± 0.'045 
1'.'582±0'.'157 
0'.'801±0'.'117 
2'.'286±1'.'411 
1'.'425±1'.'107 
2'.'761 ± 0'.'326 
0'.' 777 ±0'.' 219 
0'.'871±0'.'163 
1'.'375 ± 0'.'071 
1'.'283 ± 0'.'201 
3'.'008 ± 1'.'797 
0'.'704±0'.'126 
1'.'597 ± 0'.'057 
0'.'892 ± 0'.'087 
0'.'799 ± 0'.'026 
1'.'684 ± 0'.'081 
0'.'468±0'.'154 
1'.'904 ± 0'.'034 

1'.'360 ± ■ ■ ■ 
1'.'860 ± 0'.'260 
0'.'970 ± 0'.'382 



1.3 ±0.0 
1.5±0.1 
0.9 ±0.2 
2.0 ±0.1 

1.4 ±1.0 

1.3 ±0.1 
0.8 ±0.2 

2.0 ±0.1 
1.6 ±0.0 

1.1 ±0.1 

1.4 ±0.3 
0.9 ±0.2 
1.6 ±0.0 
1.8 ±0.2 
1.0 ±0.0 
1.6 ±0.1 
1.0 ±0.3 

1.5 ±0.0 
1.0±--- 

1.6 ±0.3 
0.8 ±0.4 



0.143 ±0.015 
0.080 ± 0.008 
0.266 ± 0.058 
0.099 ±0.136 
0.290 ±0.132 
0.348 ±0.123 
0.155 ±0.023 
0.140 ±0.084 
0.052 ± 0.059 
0.201 ± 0.004 
0.053 ± 0.004 
0.074 ± 0.002 
0.071 ±0.023 
0.280 ±0.187 
0.083 ± 0.000 
0.299 ± 0.033 
0.166 ±0.094 
0.030 ±0.016 
0.292 ±0.015 
0.087 ± 0.063 
0.021 ±0.029 



-7?9±2?8 
+64?2±3?9 
±72?4±19?7 
+4?3±16?4 
+71?7±2?5 
-28?9±0?7 
-44?9±13?0 
-30?7±32?9 
-79?3±78?5 
-54?3±6?2 
±21?5±16?0 
±80?2±--- 
-47?9±78?8 
±74?3±11?5 
-13?3±2?7 
-54?1±2?9 
-32?5±--- 
-50?2±15?6 
-10?6±19?7 
-26?2±7?1 
-51?4±43?6 



0.134 
0.171 
0.116 
0.323 
0.171 
0.142 
0.086 
0.331 
0.097 
0.117 
0.207 
0.082 
0.066 
0.361 
0.120 
0.080 
0.079 
0.107 
0.250 
0.146 
0.085 



Table 

PSF. 



6. Systematic errors due to uncertainties in the 



Quantity mean se 



median 



AWo 
ArT 

Arft 
Art 

Ae 

A^ 



0.4 

O'.'OIO 

0'.'079 

0'.'154 

0.1 

0.024 

31?8 



0.1 

0'.'003 

0'.'026 

0'.'051 

0.0 

0.004 

15? 2 



0.2 

0'.'004 

O'.'OSO 

o'.'oeo 

0.0 
0.022 

7?4 



3.1. Core Radii 

Fig. O shows the distribution of King core radii for the 
NGC 5128 GC candidates and for the Milky Way GCs. 
The Milky Way GCs have been shifted to the distance 
of NGC 5128 and their core radii have been converted 
to arcseconds (1" = 17.45 ± 0.97 pc for d = 3.6 ± 0.2 
Mpc) . Galactic GCs with a central brightness cusp (which 
is believed to be a signature of a collapsed core) have 
been omitted from the sample. The faintest GC candi- 
date in our sample has an absolute total magnitude of 
.My.tot = —6.4 ± 0.2. In order to ensure that we are com- 
paring similar objects we have excluded all Milky Way 
GCs fainter than My^tot = —6.4. 

In order to facilitate a comparison between the two 
GC systems we plotted the fraction, fi = rii/n, of the 





1 


1 




2 


<^f.^-- 


-n\ 


n. 


+ 


n 



total number of GCs in each bin, n^, where n is the total 
number of GCs in each data set. The uncertainty in Ji was 
computed from the Poisson uncertainties in the number of 
GCs in each bin, and in the total number of GCs using 



(3) 



The mean core radius for the 21 NGC 5128 GC candidates 
is f^ = 0'.'07 ± O'.'Ol (standard error) while the mean for 
the 73 selected Milky Way GCs is f^ = C'll ± O'.'Ol (se). A 
Kolmogorov-Smirnov (KS) test shows that we can reject 
the hypothesis that the two samples are drawn from the 
same distribution at the 46% confidence level. Therefore, 
there is no evidence that the core radii of the GC can- 
didates in NGC 5128 are distributed differently from the 
core radii of the Milky Way GCs. 

The most noticeable difference between the core radii 
of the NGC 5128 GC candidates and the core radii of the 
Milky Way GCs in Fig. |ll| is the lack of a tail extend- 
ing to large core radii in the NGC 5128 data. This may 
be an artifact of the small number (21) of NGC 5128 GC 
candidates in our sample. The mean of a distribution is 
sensitive to the presence of tails and outliers, but the me- 
dian is much more robust against outliers. Therefore we 
computed the median core radius for each data set. Both 
the NGC 5128 GC candidates and the Milky Way GCs 
have median core radii of [tc] = 0'.'07 (= 1.22 pc, or ^ 0.7 
pixels on the WF CCDs and 1.4 pixels on the PC CCD). 
The similarity in the median core radii suggests that the 
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Fig. 11. This figure compares the distribution of core radii 
for GC candidates in NGC 5128 with the distribution of 



core radii for selected GCs (see Sect. 'iA) in the Milky 
Way. The vertical axis is the fraction of the total number 
of GCs and the error bars show the Poisson uncertainties 
in each bin. 



"typical" core radius of a GC candidate in NGC 5128 is 
similar to that of the Milky Way GCs. 

There may be systematic biases in the core radii that 
we have derived for the GC candidates in NGC 5128. Fit- 
ting Michie-King models to objects with core radii that 
are similar to the pixel scales of the images requires that 
the centers of the objects be accurately known. Small er- 
rors in determining the center of a candidate GC, and 
small systematic errors introduced by integrating Michie- 
King model profiles over the area of a pixel, may be suf- 
ficient to bias the fitted core radii towards smaller values. 
In addition to pixelation effects, the similarity between 
the core radii and the FWHMs of the PSFs may also be 
biasing our fits toward smaller core radii. 

3.2. Tidal Radii 

The tidal radius of a GC is affected by the gravitational 
potential of its pare nt gal axy (e.g., Innanen et al 



1983 



Heggie & Ramamani 1995 ). In order to compare the tidal 
radii of GC candidates in NGC 5128 with those of GCs 
in the Milky Way it is necessary to correct for the tidal 
fields of both galaxies. The first step is to normalize the 
tidal radius of each GC candidate by its mass, Alci, to get 
Ttl M.Jy . If we assume that the gravitational potentials of 



the Milky Way and NGC 5128 can be approximated by a 
spherical logarithmic potential of the form 



<f> = y,2tln(i?2 + i?2)+in(c), 



(4) 



where R is the galactocentric distance, Rg is a scale length, 
and C is a constant, then we can compare the mean value 
of the normalized tidal radii of the GG candidates using 



rt 



M 



1/3 



Rp 



2/3 



G 



25(e) 



1/3 



(5) 



where T^ot is the amplitude of the flat part of the rota- 
tion curve of the galaxy, G is the Newtonian gravitational 
constant, and g{e) is a slowly varying function of the or- 
bital eccentricity of the GG that has values of g{0) = 1 
for circular orbits. The mean value of the normalized tidal 
radius of the GC candidates in NGC 5128, {n/MK^), can 
then be related to the mean normalized tidal radius of the 
Milky Way GCs by 



n 



- Yl 



2/3 



M 



1/3 



/ Kot,Mw \ ( gmv/{e) \ 
\ Kot ) \ 9(e) ) 



1/3 (^^2/,,^ 



{R 



J73~ 
'p,MW 



•> 



M' 



(6) 



MW 



where the subscript MW denotes the value for the Milky 
Way. 

Eq. ^ assumes that the shape of the galactic potential 
is the same in both galaxies, but allows the total mass, 
as parameterized by K-ot, of each galaxy to vary. It also 
requires a knowledge of the distribution of GG orbits in 
each galaxy, as parameterized by g{e) and Rp. We as- 
sumed a rotation velocity of T^ot.MW = 220 kms^^ for the 
Milk y Way and Kot = 245kms-i for NGC 5128 (Hui 
et al. 1995| ). The only information available on the distri- 
bution of orbits for GG candidates in NGC 5128 is the 
projected radial distances of the GG candidates from the 
center of NGC 5128, so we have assumed that the NGC 
5128 GG system is dynamically similar to the Milky Way 
GC system. This involves two assumptions about the na- 
ture of the GC orbits. First, we assume that the mean ec- 
centricity of the NGC 5128 GC orbits is the same as that 
for t he Mil ky Way GCs {e^J^ = 0.6 ± 0.1, Odenkirchen 
et al. 1997). Second, we assume that each GC is near the 
apogalacticon of its orbit, Ra, so R ^ Ra- This is a reason- 
able assumption since a GC having e ^ 0.6 spends ^ 80% 
of its time in the outer half of its orbit. An additional 
complication is that the observed galactocentric distance 
for a NGC 5128 GC candidate is the projection of the true 
galactocentric distance onto the plane of the sky. If we as- 
sume that the orbits are oriented randomly then Ra for a 
particular object's orbit is most likely to be observed to lie 
at a projected distance oi D = (8/7r^)i?a. Therefore, we 
estimated the perigalactic distances for each of the NGC 
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Fig. 12. This figure com pare s the distribution of normal- 
ized tidal radii (see Sect. 3_^) for GC candidates in NGC 
5128 with the distribution of normalized tidal radii (af- 
ter correcting for the difference in mass between the two 
galaxies) for GCs in the Milky Way. The vertical axis is 
the fraction of the total number of GCs and the error bars 
show the Poisson uncertainties in each bin. 



r^ (arcsec) 

Fig. 13. This figure compares the distribution of half-mass 
radii for GC candidates in NGC 5128 with the distribution 
of half-mass radii for a subset of GCs in the Milky Way. 
The vertical axis is the fraction of the total number of 
GCs and the error bars show the Poisson uncertainties in 
each bin. 



5128 GC candidates from the observed distance from the 
center of NGC 5128 using 



1 - e TT^ 

^ 1 + e 8 



(7) 



Fig. 02^ shows the distribution of the normalized tidal 

radii for the NGC 5128 GC candidates and the Milky Way 

GCs. The cluster masses were computed from their total 

y-band luminosities assuming a mass-to-light ratio of two. 

1 /^ 
Our sample of GC candidates has {rt/ M^[ ) ~ 0.36 ± 

0.05 (se) pc/M^(^ (iV = 21). The mean normahzed tidal 



1/3, 



/MW 



radius for 73 selected Milky Way GCs is {rt/M 
0.69±0.08 (se) pc/M^^ The multipUcative factor in Eq. | 

1/3 

is 0.703, which yields a corrected {rt/Mjy )mw of 0.49 ± 

1/3 

0.06 (se) pc/ Mq . A KS test says that we can reject the 
hypothesis that the two samples are drawn from the same 
distribution at the 74% confidence level. Therefore, there 
is insufficient evidence to state that the distribution of the 
tidal radii of the NGC 5128 GC candidates differs from 
that of the Galactic GCs if the difference in the masses of 
the two galaxies is taken into account. However, we wish to 
stress that this calculation assumes that the distribution 
of GC orbits are statistically similar for both galaxies. 



3.3. Half- Mass Radii 

The distribution of half-mass radii for the NGC 5128 GC 
candidates is shown in Fig. |l^ along with the same dis- 
tribution for our subsample of 73 Milky Way GCs. The 
half-mass radii for the Milky Way GCs were determined 
by computing a Michie-King model (with concentrations 
and core radii taken from W. Harris 1996| ) for each Milky 
Way GC. This allowed us to make a direct comparison 
between the half-mass radii of the best-fitting single-mass 
Michie-King models for the NGC 5128 GC candidates and 
the half- mass radii of the best-fitting single-mass Michie- 
King models for the Milky Way GCs. 

The mean half-mass radius for the 21 NGC 5128 GC 
candidates is f/T = 0'.'73 ± 0'.'08 (se) while the mean for 73 
selected Milky Way GCs is r^ = 0'.'67±0'.'05 (se). However, 
a KS test says that we can reject the hypothesis that the 
two samples are drawn from the same distribution at only 
the 74% confidence level, so there is no evidence that the 
distribution of half-mass radii in our sample of NGC 5128 
GC candidates is different from that of the GCs in the 
Milky Way. 



3.4- Ellipticities 

The distribution of ellipticities for the NGC 5128 GC can- 
didates is shown in Fig. |l|. The 21 NGC 5128 objects have 
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Fig. 14. The upper panel shows the distribution of ellip- 
ticities for the GC candidates in NGC 5128. The lower 
panel shows the distribution of el liptici ties for the Milky 
Way's GCs (from White & Shawl |1987D. 



e = 0.15 ± 0.02 (se) while the White & Shawl ( |1987| ) sam- 
ple of Milky Way GCs yields e = 0.07 ± 0.01 (se) for the 
73 Milky Way GCs. A KS test says that we can reject the 
hypothesis that the two samples are drawn from the same 
distribution at the 99.7% confidence level. Therefore, we 
conclude that the NGC 5128 GC candidates may have a 
different distribution of ellipticities from the Milky Way 
GCs. 

The NGC 5128 GC candidates appear to be system- 
atically more elliptical than the Milky Way GCs. There 
appears to be a lack of objects with low ellipticities and an 
excess of GC candidates with e ^ 0.3. The lack of GC can- 
didates with e < 0.05 is probably due to the elliptical PSF 
not being fully removed from the data. Another possible 
source of ellipticity is the stochastic distribution of bright 
stars near the center of the object. Geisler & Hodge (1980) 



found that the random placement of stars with respect to 
the adopted center of a GC can introduce a systematic 
error in the observed ellipticity of -1-0.045 ± 0.015 for GCs 
which are intrinsically spherical. This effect acts to make 
nearly spherical GCs appear to be more elliptical than 
they actually are. They also found that this systematic 
error decreases as the intrinsic ellipticity of the GCs in- 
creases. This would explain the lack of nearly circular GC 
candidates in NGC 5128, relative to the Milky Way. 

AM97 identified 125 GC candidates in the inner 2.'8 x 
2f8 of NGC 5128. Table @ fists the elhpticities from those 
GC candidates in common between the two studies. The 



Table 7. 


A comparison of our ellipticities witl 


AM97. 












ID 


AM97 


Our e ± CT 


AM97e 




9 


M25 


0.052 ±0.059 


0.02 




10 


M24 


0.201 ±0.004 


0.09 




12 


M23 


0.074 ±0.002 


0.04 




14 


M20 


0.280 ±0.187 


0.06 




15 


M18 


0.083 ±0.000 


0.05 




16 


M17 


0.299 ±0.033 


0.14 




18 


M15 


0.030 ±0.016 


0.09 



mean difference between our ellipticities and the AM97 
values is Ae = e(this work) - e(AM97) = -0.024 ± 0.041. 

3.5. Correlations of Cluster Properties with Distance 

GCs live within the tidal field of their parent galaxy, 
and the properties of a GC may depend on its position 
within the protogala ctic c loud at the time the GC formed 
(e.g. Murray & Lin |l992|) . In Fig. |l| we plot logio(?'c), 
log]^Q(r/,,) (both radii in arcseconds), c, and e as functions 
of \ogj^Q{D), where D is the observed distance from the 
center of NGC 5128 in arcminutes. None of the properties 
show strong correlations with the distance from the cen- 
ter of NGC 5128, but the trends are qualitative ly sim ilar 
to the trends found by Djorgovski & Meylan (1994) for 
the same properties of the Milky Way's GCs (see their 
Fig. 8). The correlations for logj^Q (re) and c for the NGC 
5128 GC candidates are in the same sense as those for 
the Milky Way GCs, although Spearman rank correlation 
coefficients suggest that the correlations are weaker for 
the NGC 5128 GC candidates. This is to be expected be- 
cause we are using the observed distance from the center 
of NGC 5128 (i.e. the distance projected onto the plane of 
the sky), while Djorgovski & Meylan (1994) used the true 
Galactocentric distance. 

We find a weak trend for the half-mass radius of a GC 
candidate to decrease as the distance from the center of 
NGC 5128 increases, while the opposite trend is seen in 
the Milky Way. The Spearman rank correlation coefficient 
for our data is —0.356, which corresponds to a significance 
of 0.887, while for the Milky Way GCs the correlation co- 
efficient is ±0.478. Therefore, there is insufficient evidence 
to for us to conclude that the observed trend for r^ to de- 
crease as the galactocentric distance increases is real. 

Djorgovski & Meylan ( |1994[ ) found no correlation be- 
tween the ellipticity of a GC and its Galactocentric dis- 
tance. We find a correlation coefficient of —0.257, corre- 
sponding to a significance of 0.739. An examination of 
Fig. na, however, suggests that there is no significant cor- 
relation between ellipticity and galactocentric distance for 
the NGC 5128 GC candidates. 



14 



S. Holland et al.: Globular clusters in NGC 5128 



-0.5 



^ -1 



O 



1.5 



-2 
0.4 
0.2 

? 



§^-0.2 



-0.4 
-0.6 



+ 0.03 (10.7%) 



I I I I I I I I 
-0.356 (88 



(a)J 

++K 

.7%); 



3 

o 2 

1 



0.4 
^ 0.2 



-0.2 0.2 
log,„(D) 



(c)- 
0.4 



-0.426 (94.6%) 



(b)- 

I I I I I I I I I I I 

-0.257 (73.9%)_ 



(d) _ 

-0.2 0.2 0.4 
log,o(D) 



Fig. 15. This figure shows the relationships between four 
properties of the NGC 5128 GC candidates and their ob- 
served distance, D, from the center of NGC 5128. The 
Spearman rank correlation coefficients, and the confidence 
levels for rejecting the hypothesis that the observed corre- 
lations are due to chance, are shown at the top of each 
panel. The lines show the best-fitting straight lines to 
the data. Panel (a) shows the core radius (in arcseconds), 
panel (b) shows the central concentration, panel (c) shows 
the half-mass radius (in arcseconds) of the best-fitting 
single-mass Michie-King model, and panel (d) shows the 
ellipticity. 

4. Colors 

4-.1. The Photometry 

Performing accurate photometry on the NGC 5128 GC 
candidates is difficult. These GC candidates are extended 
objects, not point sources, so it is not possible to deter- 
mine their magnitudes simply by fitting a scaled PSF to 
each candidate. Since each GC candidate has a unique 
size and radial profile (parameterized by the concentra- 
tion, c), aperture photometry with a single aperture will 
not return accurate magnitudes. Therefore, we elected to 
use the total flux in the best-fitting Michie-King model as 
the best estimate of its flux. 

We converted the total instrumental fluxes to the stan- 
dard Johnson-Cousins VI ma gnitudes using the prescrip- 
tion of Holtzman et al. ( 1995 ). The calibration equations 
we used are: 



V = -2.5 logio(CF555w) + (-0.052 ± 0.007)(y - I)„ 

+{0.027 ±om2)(y - 1)1 



-(21.725 ± 0.005) + 2.5 logio(G,), (8) 



/ = -2.5 logio(CF8i4w) + (-0.062 ± 0.009)(F - I)„ 
-f(0.025± 0.002) (T/-/)o 
-F(20.839± 0.006) -h 2.5 logio(G»), (9) 

where C is the count rate in ADU/second after correct- 
ing for charge transfer efficiency (CTE) effects. We ap- 



plied the CTE corrections of Whitmore & Heyer (1997) 



for a five pixel aperture. G is the gain ratio between the 
14 e~ gain state (which was used for the calibration ob- 
servations) and the 7 e~ gain state (which was used for 
the N GC 5 128 observations). As discussed in Holtzman 
et al. (1995), reddening corrections need to be apphed be- 
fore the instrumental magnitudes are calibrated to the 
standard Johnson-Cousins system. This is because the 
WFPC2 filters have different band passes and effective 
wavelengths than those of the standard V and / filters. We 
adopted a foreground reddening in the direction of NGC 
5128 of Eb-v = 0.11 ± 0.02 from the reddening maps of 
Burstein & Heiles (1982) and assumed Ey-i = 1.36£'b_v' 
(Taylor |1986| Fahlman et al. |1989D . We used the extinc- 
tion corrections for KOIII stars from Table 12 of Holtzman 
et al. ( |1995| ) to obtain Ay = 0.340 and Ai = 0.201. 

Fig. |16| shows the distribution of colors for our 21 GC 
candidates, and for 62 spectroscopically confirmed GCs 
in NGC 5128 (HG92). We converted the HG92 C - Ti 
colors toV_I colors using V-I = 0.115 + 0.514(C - Ti) 
(Geisler 1996 ). The VI magnitudes determi ned th is way 
should be treated with caution since Geisler ( |l996 ) found 
that C—Ti colors do not reproduce V—I colors particularly 
weU. There is only one object (#12, V = 20.30 ± 0.03, 
V—I = 0.40 ±0.04) which appears significantly bluer than 
the other GC candidates. Based on its color this object 
may be a young GC. Visual inspection of the WFPC2 
images (see Fig. ^) shows that this object appears to be a 
normal GC in NGC 5128. The V- and /-band photometry 
of all of our GC candidates are listed in Table g. The R- 
, J-, H-, and A'-band photometry from AM97 are also 
given. 

4.2. Reddening Within NGC 5128 

In order to estimate the amount of internal reddening 
within NGC 5128 in the direction of each GC candidate 
we determined the color of the background near each ob- 
ject. This was done by letting the mean background color 
be a free parameter during the Michie-King model fitting 
process. The expected unreddened V—I color at the loca- 
tion of each GC candidate was then subtracted from the 
observed mean background color to get an estimate of the 
internal reddening. The expected color of the background 
that each GC candidate sits on was determined as follows. 
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ID 



V 



I 



V~I 



R 



J 



H 



K 



1 


19.33 ±0.02 


18.42 ±0.02 


0.91 ±0.03 
















2 


18.21 ±0.02 


17.22 ±0.02 


0.99 ±0.02 
















3 


21.32 ±0.04 


19.94 ± 0.03 


1.38 ±0.05 
















4 


21.42 ±0.04 


19.31 ±0.03 


2.11 ±0.05 
















5 


19.76 ±0.02 


18.49 ± 0.02 


1.27 ±0.03 
















6 


19.40 ± 0.02 


17.58 ±0.02 


1.82 ±0.03 
















7 


20.91 ±0.03 


19.51 ±0.03 


1.40 ±0.04 
















8 


20.08 ± 0.03 


19.11 ±0.02 


0.97 ±0.04 
















9 


19.25 ±0.02 


18.21 ±0.02 


1.04 ±0.03 


19.17 


17.98 ±0.28 


17.62 ± 0.45 




10 


19.54 ±0.02 


18.39 ± 0.02 


1.15 ±0.03 


19.62 


17.70 ±0.19 


16.68 ±0.15 


16.56 ± 0.26 


11 


19.45 ± 0.02 


18.05 ± 0.02 


1.40 ±0.03 










12 


20.30 ± 0.03 


19.90 ± 0.03 


0.40 ±0.04 


20.51 


18.56 ±0.38 


17.74 ±0.38 




13 


18.77 ±0.02 


17.62 ± 0.02 


1.15 ±0.03 










14 


20.01 ±0.02 


19.15 ±0.02 


0.86 ±0.03 


20.05 


16.87 ±0.16 


16.62 ±0.30 


16.19 ±0.25 


15 


17.56 ± 0.02 


16.44 ± 0.02 


1.12 ±0.03 


17.39 


15.87 ±0.07 


15.13 ±0.06 


14.96 ± 0.08 


16 


18.45 ± 0.02 


17.52 ± 0.02 


0.93 ±0.03 


18.53 


17.00 ±0.12 


16.27 ±0.12 


16.66 ± 0.24 


17 


20.72 ± 0.03 


19.64 ± 0.03 


1.08 ±0.04 










18 


17.07 ±0.02 


16.08 ±0.02 


0.99 ±0.03 


17.09 


15.49 ± 0.05 


14.93 ± 0.04 


14.58 ±0.05 


19 


19.97 ±0.03 


17.78 ± 0.03 


2.19 ±0.04 










20 


20.03 ± 0.03 


17.65 ±0.03 


2.38 ±0.04 










21 


18.41 ±0.02 


17.41 ±0.02 


1.00 ±0.03 

















From Fig. 5 of van den Bergh ( |1976[ ) we derived 
B-V = 0.153 logio(i)) ± 0.924, 



(10) 



where D is the projected distance from the center of NGC 
5128 in arcminutes. We estimate that the uncertainty in 
the B — V values from Eq. 



10 



IS 



±0.1. T\ie B-V col- 
ors were converted to (V — /)q colors by subtracting the 
adopted Milky Way reddening of Eb^v = 0.11 ± 0.02 and 
using 



(t/-/)o == 0.745(B-y)o ± 0.399, 



(11) 



which we derived from the Galactic Globular Cluster Cat- 



alogue of W. Harris (1996). For each GC candidate in 



NGC 5128 we computed the expected background color 
using Eqs. 00 and O and subtracted this from the mean 
color of the unresolved background around each GC can- 
didate. The internal reddening due to dust in NGC 5128 
was then computed by subtracting the expected color from 
the observed color for each GC candidate. The internal 
reddenings, as well as the dereddened Vq and (T^ — /)q val- 
ues for each GC candidate, are listed in Table y. The Vb 
and {V — I)(y values are corrected for both the Burstein & 
Heiles ( |1982 ) reddening and our estimate of the internal 
reddening in NGC 5128. The uncertainty in each redden- 
ing estimate is ^ 0.08 and negative internal reddenings 
were set to Ey-i = 0. 

In order to check these differential reddenings, we re- 
peated our calculation using the color of a typical galaxy 
with the same morphological type as NGC 5128. The 



Table 9. The estimated internal reddening within NGC 
5128 along the line of sight, the dereddened V^-band mag- 
nitude, and the dereddened (F — /)q color for each GC 
candidate. 



ID 


Ev-i 


Vo ± avo 


(1/-J),±a(,._,)„ 


1 


0.17 


18.94 ±0.06 


0.74 ± 0.08 


2 


0.15 


17.87 ±0.06 


0.84 ± 0.08 


3 


0.06 


21.18 ±0.07 


1.31 ±0.09 


4 


0.00 


21.42 ±0.07 


2.12 ±0.10 


5 


0.00 


19.76 ± 0.06 


1.27 ±0.09 


6 


0.07 


19.25 ± 0.06 


1.75 ±0.09 


7 


0.14 


20.59 ± 0.07 


1.27 ±0.09 


8 


0.06 


19.94 ± 0.06 


0.90 ± 0.09 


9 


0.08 


19.08 ± 0.06 


0.96 ± 0.08 


10 


0.11 


19.28 ± 0.06 


1.03 ± 0.09 


11 


0.09 


19.26 ± 0.06 


1.32 ± 0.09 


12 


0.14 


19.98 ± 0.06 


0.26 ± 0.09 


13 


0.10 


18.53 ± 0.06 


1.04 ±0.08 


14 


0.10 


19.77 ± 0.06 


0.76 ± 0.09 


15 


0.06 


17.42 ± 0.06 


1.06 ± 0.08 


16 


0.14 


18.13 ±0.06 


0.79 ± 0.08 


17 


0.13 


20.42 ± 0.07 


0.95 ± 0.09 


18 


0.10 


16.83 ± 0.06 


0.89 ± 0.08 


19 


0.10 


19.76 ±0.07 


2.10 ±0.09 


20 


0.00 


20.03 ± 0.07 


2.38 ± 0.09 


21 


0.22 


17.91 ±0.06 


0.78 ± 0.08 



Third Refere nce C atalogue of Bright Galaxies (de Vau- 
couleurs et al 1991) lists NGC 5128 as being an intermedi- 
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Fig. 16. The upper panel shows the distribution of colors 
for our 21 GC candidates in NGC 5128 while the lower 
panel shows the distribution of colors for 62 spectroscopi- 
cally confirmed GCs in NGC 5128 from HG92. The colors 
have been corrected for reddening in the Milky Way in the 
direction of NGC 5128, but not for internal reddening in 
NGC 5128 itself. 



ate SO galaxy with a morphological type of T = —2 ± 0.3. 
Buta & Williams ( 1995| ) find a mean color for T = —2 
galaxies of (V^ — /}totai = 1.145 ± 0.069 based on a sample 
of 55 galaxies. This is within Itr of the expected back- 
ground colors derived from the observed color gradient in 
the outer regions of NGC 5128. If we assume that the un- 
resolved light from NGC 5128 is the same as that from a 
typical T = —2 galaxy, but has been reddened due to the 
presence of dust from the merger, then the excess redden- 
ing can be estimated by subtracting the mean color of a 
T = — 2 galaxy from the observed color of the unresolved 
background around each GC candidate. The differential 
reddenings that we obtain by assuming that NGC 5128 is 
a T = —2 galaxy are consistent with those obtained using 
the van den Bergh ( 1976 ) color gradient. 

The first problem with this method of estimating the 
internal reddening in NGC 5128 is that NGC 5128 is not a 
normal elliptical galaxy, but an elliptical galaxy that has 
undergone a merger with a small late-type spiral galaxy. 
Therefore, the unresolved light near the center will be a 
combination of light from the original galaxy and from 
stars in the captured spiral galaxy. In our derivations of 
the internal reddening we have assumed that the only 
changes in the color of the central regions of NGC 5128 
are those due to dust, and ignored changes in color due 
to star formation and differences in the underlying stellar 



population. Buta & Williams (1995) list mean colors for 
late spiral galaxies of (F — /)totai ~ 0.8 to 0.7, depending 
on the morphological type of the spiral. Therefore, the 
contribution from the merged spiral galaxy will cause the 
actual unreddened color of the unresolved background in 
NGC 5128 to be as much as a few tenths of a magnitude 
bluer than we have assumed. This will lead us to underes- 
timate the internal reddening within NGC 5128 by a few 
tenths of a magnitude. It is possible that some of the GC 
candidates are being seen against regions of recent star 
formation in NGC 5128. This would also lead us to un- 
derestimate the internal reddening by a few tenths of a 
magnitude. 

A second problem is that the structure of the dust 
features in NGC 5128 changes on spatial scales of less 
than ^ 1". The mean tidal radius of the GC candidates 
is r7 = 1'.'42 ± 0'.'15 (se), and the surface brightness of 
the background was determined while fitting the Michie- 
King models by taking the mean value of the pixels that 
fell within the fitting box (64 x 64 pixels) but beyond 
the tidal radius. Therefore, small-scale spatial structure 
in the dust lane could introduce errors of several tenths of 
a magnitude in our estimates of the internal reddening in 
NGC 5128. 

Finally, the reddening corrections are not being made 
in the WFPC2 photometric system. This may introduce 
small errors in the reddening corrections that we deter- 
mine. We believe, however, that these errors will be small 
compared to other uncertainties in the method. In light 
of these uncertainties we believe that it is possible that 
our estimates of the differential reddening for each GC 
candidate are only accurate to ~ ±0.3 mag. 

Fig. |7| shows the color distribution of the GC candi- 
dates after correcting for internal reddening within NGC 
5128. There appear to be three populations in the up- 
per panel of Fig. 17, The largest population is centered 
at (F — /)q ^ 1.0 and contains 16 of the 21 GC candi- 
dates. This population appears to be similar to the spec- 
troscopically confirmed GCs of HG92 and probably rep- 
resents a genuine old population of GCs in NGC 5128. 
The similarity between the colors of these GC candidates 
and the dereddened colors of the Milky Way GCs suggests 
that this population is not heavily reddened and therefore 
probably lies on the near side of NGC 5128. The second 
component is the four objects with {V — I)^ > 1.4. These 
objects are probably GCs that are being seen through 
large amounts of dust in NGC 5128. Alternately, they 
could be background galaxies that have been misidenti- 
fied as GC candidates. The third component is the single 
blue object with {V — I)^ = 0.26 ± 0.09. An examination 
of the image of this object (Fig. ^ suggests that it is a 
legitimate GC candidate so its blue color makes it the best 
candidate for being a young GC in our sample. 
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Fig. 17. The upper panel shows the distribution of col- 
ors for the GC candidates in NGC 5128 after correcting 
for the estimated internal reddening due to dust in NGC 
5128. The middle panel shows the distribution of colors 
for the GC candidates without any correction for inter- 
nal reddening within NGC 5128 (but with a correction for 
Galactic reddening in the direction of NGC 5128). The 
lower panel shows the distribution of dereddened colors 
for the Milky Way GC system. 

4-. 3. Iron Abundance 

The iron abundance for each GC candidate was estimated 
using 



[Fe/H] 



-6.568 + 6.733(y-/)o - 1.372{V-I)l, (12) 



which was determined from 45 Galactic GCs with low 
reddenings and good metallicity estimates, and from 13 
metal-rich GCs in NGC 1399 (Kissler-Patig et al. 



1998) 



If the blue color of object #12 is due solely to its 
iron abundance, then it has [Fe/H] — —4.91 ± 0.54. 
This implausibly low iron abundance suggests that age 
is primarily responsible for at least some of the blue 
color, supporting the notion that #12 is a young ob- 
ject. The mean iron abundance of the 16 GC candidates 
with 0.6 < {V~I)q < 1.4 is [Fe/H] = -1.27 ± 0.20 (sc), 
which is within 2a of the [Fe/H] = —0.8 ± 0.2 value found 
by HG92. This suggests that the colors of these objects 
are consistent with them being old GCs, or moderately 
reddened intermediate-age GCs. The four objects with 
{V — I)q > 1.4 all have iron abundances of greater than 
the Solar value ([Fe/H] = +1.44 ± 0.15). This suggests 
that the red colors of these objects are primarily due to 
dust along the line of sight within NGC 5128, not high 



iron abundances, although we can not rule out the possi- 
bility that at least some of these GC candidates have high 
iron abundances. 



5. Young GCs in NGC 5128 

We would expect young GCs to be brighter and bluer than 
their old counterparts, although reddening uncertainties 
complicate the interpretation of the observed colors for the 
NGC 5128 GC candidates. The mean {V-I)^ color of the 
97 Milky Way GCs that have both V-I and Er-v values 



listed in the catalogue of W. Harris ( |1996D is 0.91 ± 0.09 
(standard deviation). This is consistent with the colors of 
the majority of the GC candidates that we found in the 
inner regions of NGC 5128. 



Schweizer & Seitzer (1993) find that young GCs in 
NGC 7252 have V — I colors that are ~ 0.5 mag bluer 
than old GCs, which suggests that objects in NGC 5128 
with colors bluer than V~I ^ 0.8 may be young (ig < 1 
Gyr) GCs. We find five objects with (V-I)^ < 0.8. One 
of these, object #12, is unusually blue, with {V — I)q = 
0.26 ± 0.09. AM97 identified this object as M23 and mea- 
sured a J — H color of 0.83. A visual examination of the 
HST images (see Fig. 0) shows nothing unusual about this 
object. If we assume that the blue color of this object is 
due to it being a young star cluster, we can use the models 
of Bruzual & Chariot ( 1993 ) to estimate its age and mass 
from its {V — I)q color and its integrated V magnitude. 
Fig. [l8fi shows the evolution of the (¥ — 1)^ color of star 
clusters for several iron abundances. This figure suggests 
an age for #12 of less than ~ 100 Myr, while Fig. |l|b 
suggests that, if it is a young, metal-poor GC, it will fade 
by between ~ 4 and ~ 6.5 mag over the next 12 Gyr. 
Therefore, when object #12 is as old as a typical Milky 

Way GC (~ 12 Gyr), it will have My between 3.5 

and —1, corresponding to a mass of between ^ 4000 and 
^ 400 Solar masses, assuming a mass-to-light ratio of two. 
W. Harris (1996) lists four Galactic GCs with masses less 
than 4000 Solar masses: AM 1 with a mass of 700 Solar 
masses. Pal 1 with a mass of 1000 Solar masses, E 3 with 
a mass of 2100 Solar masses, and Terzan 1 with a mass 
of 3400 Solar masses, although these mass determinations 
are highly uncertain. In light of the small mass for object 
#12 it is also possible that this object is a massive open 
cluster. 

The sixteen GC candidates with 0.7 < (T^-/)o < 1-4 
have colors that are consistent with being young (< 1 
Gyr), intermediate-aged (1-4 Gyr), or old (> 4 Gyr) GCs 
depending on the iron abundances and reddenings that 
are assumed. If these GC candidates have iron abundances 
similar to those found in the Galactic GCs ([Fe/H] < 0) 
then Fig. p^ suggests that their observed (y — /)q colors 
are consistent with them having ages of > 1 Gyr. If this is 
the case then Fig. Illb suggests that these objects will fade 
by between ^ and 2 mag in V by the time they are ~ 12 
Gyr old. This corresponds to these GC candidates having 
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The red colors, and high implied iron abundances, of 
the four GC candidates with {V — T)^ > 1.4 suggest that 
these objects are heavily reddened. If we assume that all 
of these objects are young, metal rich clusters that formed 
at the time of the merger event then differential redden- 
ings of Ey-i ^ 1.1 to 1.8, in addition to those derived 



Fig. 18. The upper panel shows the time evolution of 
{V — I)q for GCs with a Salpeter Initial Mass Function 
for a range of cluster masses from 0.1 to 125 Solar masses, 
and a range of iron abundances from [Fe/H] = —1.64 to 
[Fe/ H] = +1.01, using the models of Bruzual & Char- 
lot ( 1993 ). The heavy vertical lines are at 160 Myr and 500 
Myr respectively (the interval that the most recent merger 
most likely occurred in). The colors of the 21 GC candi- 
dates are plotted as squares (Galactic reddening removed) 
and circles (Galactic and internal reddening within NGC 
5128 removed). The lower panel shows the time evolution 
of the V magnitude per Solar mass using the same models. 



masses between ~ 40 000 and ~ 4 000 000 Solar masses, 
which is consistent with the masses of Galactic GCs. If, on 
the other hand, the objects have iron abundances greater 
than the Solar value then their colors are consistent with 
their being young GCs. In this case Fig. ^8|b suggests that 
they will fade by ~ 3 to 4 mag in V over then next 12 
Gyr, which would give them masses of between ^ 5000 
and ~ 250 000 Solar masses. 

These results assume that the differential reddenings 
derived in Sect. ^ are uncertain by ~ ±0.3 mag, and 
may contain systematic uncertainties of a few tenths of a 
magnitude. The {V — I)^^ colors of the sixteen GC candi- 
dates with 0.7 < (y — /)g < 1.4 are most consistent with 
young GCs if the differential reddenings have been sys- 
tematically underestimated by between 0.1 and 0.7 mag. 
If our estimates of the differential reddening for these ob- 
jects are correct then their {V — I)^ colors suggest that 
these GC candidates are normal, old, metal-poor GCs. 
Spectroscopic observations are needed to break the de- 
generacy between reddening and metallicity to accurately 
determine the ages of these objects. 



in Sect. .^.1, are required to account for their observed 
red colors. This corresponds to between ^ 2.4 and ~ 3.9 
mag of extra extinction in the V band. If, on the other 
hand, these objects are young metal poor clusters then 
the amount of additional differential reddening required 
to account for their observed red colors is between '^1.3 
and '^2.1 mag. It is unlikely that we have underestimated 
the differential reddening by this amount. If we assume 
that these objects are old GCs (~ 12 Gyr) seen through 
dust then we have underestimated the amount of differ- 
ential reddening by between r^ 0.4 and ~ 1.4 mag. It is 
possible that a combination of random and systematic un- 
certainties in the differential reddening calculations could 
result in us underestimating the reddenings of these ob- 
jects by several tenths of a magnitude, so we can not rule 
out the possibility that these four objects are old GCs. An- 
other possibility is that they are background red elliptical, 
or dwarf elliptical galaxies with concentrations similar to 
those of the Galactic GCs. Spectroscopic observations are 
needed to confirm the nature of these objects. 

In summary, the dereddened (T^ — /)q color of object 
#12 is consistent with it being either a small young GC 
or a massive open cluster. The sixteen objects with 0.8 < 
{V — I)q < 1.4 have colors that are consistent with them 
being either young GCs which suffer from up to ~ 0.7 
mag of differential reddening or old GCs that suffer from 
up to ^ 0.2 mag of differential reddening. Spectroscopic 
observations will be needed to determine the ages of these 
GC candidates. The four objects with {V — I)q > 1.4 are 
either old heavily reddened GCs or background galaxies 
that have similar concentrations to Galactic GCs. 



6. Conclusions 

We have used the HST /WFPC2 to identify 21 GC candi- 
dates in a 25n' area centered on the nucleus of the near- 
est giant elliptical galaxy, NGC 5128. There is strong ev- 
idence that this galaxy has undergone a merger with a 
small late-type spiral between ~ 160 and 360 Myr ago. We 
fit two-dimensional, PSF-convolved, single-mass Michie- 
King models to each GC candidate and derived core radii, 
tidal radii, half-mass radii, ellipticities, and position angles 
for each object. We assumed that the GC candidates were 
structurally similar to those in the Local Group; i.e., they 
could be well fit by Michie-King models, and had similar 
half-mass radii and ellipticities. Therefore, only objects 
with structural parameters similar to those of GCs in the 
Local Group were accepted as being GC candidates. It is 
possible that we have missed NGC 5128 GC candidates 
that are significantly less centrally concentrated than the 
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GCs in the Local Group. However, since we were primar- 
ily interested in using the colors of the GC candidates to 
estimated their ages we preferred to risk rejecting legiti- 
mate GC candidates rather than risk having our sample 
contaminated with stars and background galaxies. Within 
this constraint we find no evidence that the NGC 5128 GC 
candidates have core-, half-mass-, or tidal radii that are 
distributed differently from those GCs in the Milky Way 
that do not exhibit central brightness cusps. There is weak 
evidence that NGC 5128 GC candidates are systematically 
more elliptical than are the Galactic GCs. 

We have obtained V- and /-band photometry for all 
21 of our GC candidates. We find no evidence for the 
bimodal color distribution observed among GCs at larger 
distances from the center of NGC 5128 (HG92; Zepf & 
Ashman 1993), although this is likely due to the small 



sample size (21 objects), the poor metallicity sensitivity 
of the V — I color index, and confusion due to differential 
reddening within NGC 5128. 

We have identified one very blue GC candidate (#12) 
with {V — I)q = 0.26 ± 0.09. Using a color-age relation 
derived from the models of Bruzual & Chariot (1993), we 
estimate that this object has an age of less than ~ 100 Myr 
and a mass of between ~ 400 and ~ 4000 Solar masses, 
which is barely consistent with this object being a small 
GC that formed during the merger event. There are six- 
teen objects with 0.7 < {V — I)^ < 1.4 that have colors and 
integrated magnitudes that are consistent with their be- 
ing either young GCs that formed during the merger event 
or old GCs similar to those found in the Milky Way, de- 
pending on what we assume about their iron abundances 
and differential reddenings. The colors of these GCs are 
very similar to those measured by HG92, but the amount 
of reddening for each object is uncertain by ~ 0.3 mag. 
There are four GC candidates with {V — I)q > 1.4, which 
imply either implausibly high metallicities or very large 
difi^erential reddenings. We can not rule out the possibil- 
ity that some of these four objects are background galaxies 
seen through the central regions of NGC 5128. 

Spectroscopic studies will be needed to determine un- 
ambiguously if any of the central GC candidates in NGC 
5128 are young GCs that may have formed as a result 
of the merger. If young GCs in NGC 5128 can be unam- 
biguously identified spectroscopically, and their ages de- 
termined, this will make it possible to estimate the amount 
of time required for GCs to form after a galactic merger 
has occurred. 
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